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Overview 



• NASA’s space radiation transport code, HZETRN, has been used 
extensively in a variety of applications 

- Research 

- Vehicle design and optimization 

- Operations and risk assessment 

• HZETRN is a deterministic model allowing a converging sequence of 
physical approximations to be implemented 

- Direct solution of the Boltzmann equation using method of characteristics 

- Straight-ahead approximation (extensive verification and validation) 

- Bi-directional transport extension (extensive verification, limited validation) 

- 3D transport extension (limited verification) -> present focus 

• Transport extensions beyond straight-ahead approximation based on 
forward/isotropic interaction model 

- Neutron production terms can be represented as a sum of forward and isotropic 
components 

- 3D extension allows for integration of full angular dependence (future work) 
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3DHZETRN Formalism 


• Neutron production cross section is separated into forward and isotropic 
components 

• Forward component of flux is solved within the straight-ahead 
approximation 

- at any point within the geometry 

• Forward flux generates isotropic neutron source term 

- at any point within the geometry 

• Isotropic neutron field described with 
N stream directions 

- Bi-directional neutron transport (N=2) 
implemented along anti-parallel rays 

- Final step evaluates/transports light ion 
target fragments produced from isotropic 
neutron field 
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3DHZETRN Formalism 




Plot shows neutron production cross 
section for 500 MeV protons on 
aluminum 

- Forward (for) component extends up to 
high energies 

- Isotropic (iso) component associated with 
lower energy knockout and evaporation 
processes 



Plot shows finite geometry with 
boundary condition directed along Q 0 

- Forward solution at target point and at 
points along stream direction, Q, obtained 
within straight ahead approximation 

- Discrete number of Q directions (streams) 
chosen to represent isotropic field 


fin 
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Verification: Sphere benchmark (I) 




First benchmark is for a sphere of aluminum (radius 20 g/cm 2 ) with external 
boundary condition applied uniformly down onto sphere along z-axis 

- 1 956 Webber SPE spectrum 

- Solar minimum GCR ions: Z = 1, 2, 6, 26 


Convergence testing performed with N (number of streams) 

- Reasonable convergence met after N = 1 8 


Boundary condition 

^ ^ ^ 

+z 



40 g/cm 2 



Benchmark: Webber SPE 


• Fluctuations in MC results associated with ENDF elastic resonances 

• 3DHZETRN (A/=18) significant improvement over straight-ahead 
approximation (A/=1) 



Boundary condition 
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1956 Webber SPE boundary condition 



Benchmark: GCR hydrogen 



Boundary condition 


• Good agreement amongst codes for secondary neutrons 
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1977 solar min GCR hydrogen boundary condition 
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Benchmark: GCR iron 



Good agreement amongst codes for secondary nucleons 

Some disagreement in secondary alphas - driven by nuclear cross 
section models, not transport approximations 
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Verification: Sphere benchmark (II) 



• Next benchmark is for a sphere of ICRU tissue (radius 15 g/cm 2 ) surrounded by 
20 g/cm 2 thick aluminum shield with external boundary condition applied 
uniformly down onto sphere along z-axis 

- 1 956 Webber SPE spectrum 

- Solar minimum GCR ions: Z = 1, 2, 6, 26 


• Convergence testing performed with N (number of streams) 

- Reasonable convergence met after N = 22 


Boundary condition 
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Benchmark: Webber SPE 


• Straight-ahead approximation has incorrect neutron 
spectral shape 

• 3DHZETRN (N=22) in close agreement with MC 
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Benchmark: GCR hydrogen 



• 3DHZETRN (N=22) in close agreement with MC 
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Benchmark: GCR iron 



• 3DHZETRN (N=22) in close agreement with MC for 
secondary nucleons 

• Differences in secondary alphas driven by nuclear cross 
section models 
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Summary and Future Work 

• Newly developed 3DHZETRN code includes 3D corrections for 
neutrons and neutron-induced light ions 

• Benchmark comparisons with MC codes showed significant 
improvement in 3DHZETRN results and close agreement for 
nucleons 

- Non-trivial differences in nuclear physics models still present 

• Current code can work from limited set of combinatorial geometries 
(currently extending) 

• Extend to engineering ray trace representation of complex 
geometries 

- Connects to current methodologies used for vehicle design and shield 
optimization 

• Extend to 3D light ion/neutron using isotropic/straight-ahead 
interaction 

• Begin space flight validation 
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Backup 



Convergence Testing (I) 



• Image at right shows the ray-distributions 
used to evaluate isotropic particle field 

• Plot below shows nucleon and alpha 
fluences induced by Webber SPE in 
benchmark geometry 

• ISM 8 and N=22 are clearly converged 






1956 Webber SPE boundary condition 
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Benchmark: GCR alpha (I) 



Good agreement amongst codes for primary alphas and secondary 
nucleons 
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1977 solar min GCR alpha boundary condition 


Benchmark: GCR carbon 


Good agreement amongst codes for secondary nucleons 

Some disagreement in secondary alphas - driven by nuclear cross 
section models, not transport approximations 
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Convergence Testing (II) 



• Image at right shows the ray-distributions 
used to evaluate isotropic particle field 

• Plot below shows nucleon and alpha 
fluences induced by Webber SPE in 
benchmark geometry 

• ISM 8 and N=22 are clearly converged 

• ISM gives incorrect spectral shape 





20 




• 3DHZETRN (N=22) in close agreement with MC for 
secondary nucleons and primary alphas 
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Benchmark: GCR carbon 


• 3DHZETRN (N=22) in close agreement with MC for 
secondary nucleons 

• Differences in secondary alphas driven by nuclear cross 
section models 
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